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Abstract. - We propose to entangle macroscopic atomic ensembles in cavity using EPR- 
correlated beams. We show how the field entanglement can be almost perfectly mapped onto 
the long-lived atomic spins associated with the ground states of the ensembles, and how it can 
be retrieved in the fields exiting the cavities after a variable storage time. Such a continuous 
variable quantum memory is of interest for manipulating entanglement in quantum networks. 



Entanglement is one of the most intriguing feature of quantum mechanics and, since the 
enunciation of the famous Einstein-Podolsky- Rosen paradox [1], has always attracted a lot 
of attention. In particular, it is at the heart of quantum communication and quantum in- 
formation protocols such as quantum cryptography, teleportation, dense coding, quantum 
computing [2]. The past few years have seen many realizations of entangled beams in the 
continuous variable regime, using process in optical parametric amplifiers (OPAs) [3-6], 
Kerr effect in optical fibers [7,8] or in cold atoms [9]. Efficient sources of entangled beams 
now exist and strong correlations have been achieved over rather broad bandwidth [6]. In 
order to build quantum communication networks in which light beams connect atomic ensem- 
bles, a major issue is to be able to store entanglement into the atoms [10,11]. Entanglement 
between two atomic ensembles has been successfully demonstrated by Julsgaard et al. by 
sending pulses of coherent light through two atomic vapor cells [12] and measuring the outgo- 
ing field. However, the possibility to store entanglement between quantum-correlated beams 
into atoms remains to be demonstrated. In this Letter we propose a cw scheme to achieve 
entanglement between two cold atom ensembles placed in cavities by using EPR-cntangled 
beams, as produced by OPAs for instance, and coherent control fields. The entanglement 
between the beams is mapped onto the ground state spins of the atoms and no measurement 
of the field is required. Given the long lifetime of the cold atoms spin the entanglement can 
thus be stored for a rather long time when the control field is switched off. It can then be 
retrieved in the vacuum modes exiting the cavities by switching on the control field again 
after a variable storage time. We then give a method to directly measure the entanglement 
of the outgoing beams - and, consequently, the atomic entanglement - in one simultaneous 
measurement of the EPR variances with two homodyne detections and a single local oscillator. 
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Fig. 1 - (a) Entanglement mapping onto the atoms: the incoming vacuum fields A™ 2 are entangled, 
Oi ; 2 are intracavity coherent control fields. The insert shows the A structure of the atomic levels con- 
sidered, (b) Readout: the incident vacuum fields are in a coherent vacuum state, the outgoing vacuum 
fields are entangled when the control fields are switched on again. PBS: polarizing beamsplitter, LO: 
local oscillator, PZT: piezo-electric ceramic. 



As shown in fig. we consider two identical sets of ./V A-type 3-level atoms, each set 
interacting with a control field f2j, (i = 1,2) and with one of the EPR-entangled beams A{. 
Without loss of generality we assume the control fields to be er+-polarized and the entangled 
vacuum fields er_-polarized. The entangled vacuum fields can be obtained for instance from 
OPAs [5] or from an OPO below threshold [6], although other schemes can be equivalently 
envisaged. To simplify, the entanglement bandwidth is assumed to be larger than the cavity 
bandwidth re. First, we examine how to create entanglement between the atoms and start by 
studying the interaction of light with one ensemble. In previous works [13-15], we showed 
that "Raman" - or "EIT"-type interaction of light with A-type atoms could lead to squeeze the 
atomic spin, cither in a non linear regime, when the incoming field is in a coherent state ("self 
spin squeezing"), or by transfer, when the incoming field is a broadband squeezed vacuum. 
In the latter, we examined how squeezing could be transferred from the field to the atoms 
in different configurations. A quasiperfect transfer is predicted either in an "EIT" (on one- 
and two-photon resonance) or a "Raman" (on two-photon resonance, but large one-photon 
detunings) configurations. Basing ourselves on these results we consider for instance an EIT 
situation in which the fields are both one- and two-photon resonant (A = 5 = 0) and cancel 
the cavity detuning. Since (Af l ) = 0, the atoms are pumped by control field ill into level 
|2) and the spin is aligned along z in steady state: (J z i) = N/2. The collective atomic spin 
can then be treated as a harmonic oscillator, the non-commutating spin components J xl and 
J y i in the plane orthogonal to the mean spin playing the same role as the field quadrature 
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operators X\ = A\ + A[ and Y\ = i{A[ — A\). We thus seek to map the EPR correlations 
existing between the quadratures Xj, onto the two ensembles spin components J X i, J y i. 
It is possible to choose the control field pumping rate Te = ^1/7 so that the ground state 
observables evolve slowly with respect to the fields or the optical dipoles and get simplified 
equations for the atomic spin fluctuations in the Fourier domain [15] 

(Tb-tw^iH = -/3SXi n (uj) + f xl (1) 
(%-iu)5J vl (u) = -fiSYpiu) + f vl (2) 

where 70 = 70 + r^/(l + 2C) is the effective atomic decay constant, satisfying 70 -C 70 <C 7, K, 
C = g 2 N/T~f is the standard cooperativity parameter, T the coupling mirror transmission, g 
the atom- field coupling constant and f x \, f v \ are Langevin operators, the correlation functions 
of which can be calculated via the quantum regression theorem. These operators account for 
the noise due to loss of coherence in the ground state (oc 70) and for the noise contribution 
of the optical dipole via spontaneous emission (oc Te)- P = gNCli/jy/T(l + 2C) represents 
the effective coupling with the incoming EPR field A™ . To simplify we assume a symmetrical 
configuration for ensemble 2 (Cl 2 = f2i, same number of atoms, (J z2 ) — N/2, etc), so that 
the equations for the ground state spin fluctuations are similar to eqs. lllll'l) by substituting 
subscript 1 by 2. We then obtain very simple equations for the fluctuations of the ground 
state spins operators 

(7 -iu)(SJ xl -5J x2 ) = -p(6Xr-SX l 2 n ) + f xl -f x2 (3) 
(7o-ioj)(5J y i+SJ y2 ) = -[3(6Yr + 5Y 2 m ) + f y i + fy2 (4) 

which are valid if the effective optical pumping rate satisfy 70 <C 70 <C 7. Since the fluctu- 
ations of both spins are related to those of the incident beams quadratures, one expects that 
the field correlations will reflect on the atoms. 

To quantify the entanglement we make use of the inseparability criterion derived by Duan 
et al. and Simon [16], which rely on the sum of the variances of EPR-type Gaussian operators, 
such as Xx — X 2 and Y\ + Y 2 for the field. 

I f = ^[A 2 (X 1 -X 2 ) + A 2 (Y 1 +Y 2 )} <2 (5) 

is a necessary condition for modes 1 and 2 to be entangled. Moreover, when modes 1 and 
2 are symmetric it has been shown by Giedke et al. that this condition is also sufficient 
and directly related to the Entanglement of Formation (EoF), thus providing a good measure 
of entanglement [17]. For OPAs pumped below threshold strong amplitude correlations and 
phase anti-correlations exist, and both A 2 (Xi — X 2 ) and A 2 (Yi +Y 2 ) can be strongly reduced 
below the separable beams value of 2 on a broad bandwidth [6]. In a similar fashion the 
EPR-type operators for the atomic ensembles are J x \ — J x2 and J y \ + J y2 , which satisfy 
([Jxi — Jx 2 , Jyi + Jy 2 ]) = i(Jzi — Jz 2 ) = 0. The atomic ensembles are then entangled if the 
following inseparability criterion is satisfied [12] 

A 2 (J xl -J x2 ) + A 2 (J yl + J y2 ) < \(J z i)\ + \(J z2 )\ = N (6) 

It is then obvious from eqs. il-SI II) that any entanglement between modes 1 and 2 will be 
transferred to the atomic spin, provided the coupling (3 is large enough with respect to the 
noises of the process / Q j. This is indeed possible in EIT or Raman situations for which the 
noises are substantially reduced by the cooperative behavior of the atoms [15], while the 
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Fig. 2 - (a) Left axis: Inseparability criterion for the incident fields If [plain] and the atoms T a t 
[dotted: simple model approximation, dashed: full calculations] vs EPR correlations 2— Tf. Right axis: 
entanglement of formation for incident fields and atoms vs EPR correlations [parameters: C = 100, 
K = 27, 7 = IOOO70, T_E = 157]. (b) Optimized mapping fidelity rj — f(T a t)/ f(Tf) vs cooperativity. 
For each value of C the pumping rate is optimized [7 = IOOO70, k = 27,1/ = 1]. 



coupling is enhanced by the number of atoms, as we will show further. In other words, the 
signal-to-noise ratio of the transfer process is increased by the atoms cooperative behavior, 
even though the strong coupling regime is not reached [18]. Assuming broad bandwidth 
amplitude correlations and phase anti-correlations and in the regime 70 <C 70 7, k, the sum 
of the atomic EPR variances, T at = jf[A 2 (Jxi — Jxi) + A 2 (J y i + J y2 )] (normalized to 2) is 
directly related to the amount of EPR-type correlations of the incident beams Tf 
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The first term in 0, proportional to Tf, can be understood as the atom- field coupling factor 
and can be very close to 1 for C and T e large enough. The second and third terms represent 
the noises of the transfer process due to the loss of coherence in the ground state (oc 70) and 
spontaneous emission (oc Te), respectively. Both can be made small in the regime chosen for 
the pumping and for large values of C. Typical experimental values of C of 100-1000 ensure a 
quasiperfect entanglement mapping from fields to atoms. Note also that, if the incident fields 
are not entangled (Tf — 2), so are the atoms (T at — 2). The accuracy of our simple picture 
has been checked by full calculations involving the three-level atomic structure and the exact 
covariance matrix of the whole atom-field system. The results are represented in fig. [2Ja) 
in which we plot T at as a function of the EPR correlations 2 — If. Very good agreement is 
found between the simple model [Eq. 0, dotted] and the numerical simulations [dashed]. As 
a comparison the inseparability criterion for the incident fields Tf is also represented [plain] . 
In a symmetrical configuration the inseparability criterion Q is a measure of entanglement 
via the entanglement of formation [17]. The incident EPR fields and the atomic EoFs (see [19] 
for details) are also both plotted in fig. El[a) . To quantify the efficiency of the mapping we 
plot in fig. [2Kb) the ratio of the atomic EoF to that of the EPR fields, 77 = / (T at ) / f (I /) , 
as a function of the cooperativity. An excellent mapping (77 ~ 100%) is obtained for easily 
accessible values of C. 



Once the field entanglement has been mapped onto the atoms the fields can be switched 
off and the atoms return to the coherent spin state very slowly, on a time scale given by 1/70- 
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We now address the readout problem and show that it is possible to retrieve the entanglement 
in the vacuum fields exiting the cavities by switching on again only the control fields. Indeed, 
starting now with correlated atomic ensembles and incoming coherent vacuum fields we ex- 
pect the fluctuations of J x i — J X 2 and J y \ + J y2 to imprint on the fluctuations of, respectively, 
■g out _ x° ut and Y° ut + Y 2 out . As represented in fig. ^the outgoing vacuum fields are combined 
using quarter-wave plates on the first PBS while the outgoing control fields are discarded. Af- 
ter PBS1 one disposes of two entangled vacuum modes with orthogonal polarization. The 
path dephasing can easily be cancelled using the control fields interference signal to lock the 
piezo-electric ceramic. Using the method developed in ref. [20] we then rotate the polarization 
basis so as to retrieve two squeezed modes for the same quadrature. This can be easily done 
with a quarter- wave plate at 0° in order to rotate the noise ellipsoid of one mode by ir/2 with 
respect to the other. One then uses a half-wave plate at 22.5° to get the ±45° polarization 
modes, which are now squeezed for the same quadratures. To simultaneously measure the 
squeezing of both modes we use the technique of ref. [9] and combine on PBS2 the beam to be 
measured with a local oscillator polarized at 45° to the cube axes. We then perform two bal- 
anced homodyne detections, the first of which measures the noise of X[ = (X° ut — X 2 ut ) / \[2 1 
the second, X' 2 = {Y° ut + Y 2 out )/y/2. The sum of the signals therefore gives directly the value 
of the inseparability criterion JSJ. 



More quantitatively one may express the correlation function of X[ in the regime 70 <C 
70 < k, 7 as functions of the atomic EPR variances at the switching time 



C[{t,t') = (SX[(t)SX[(t')) = 6(t - t') 



ACT, 



(i + 2cy 



1 - 
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A similar expression holds for the correlation function of X 2 , replacing J x \ — J x2 by J y \ + J y2 . 
When the atoms are in coherent spin states (A 2 (Jj;i — J x a) = N/2) one retrieves the standard 
(5-function of a free field. We must now specify how the homodyne detections are performed. 
In order to correctly measure the squeezing of modes A[ and A' 2 , and, therefore, the field 



entanglement at the output of the cavity, we choose for the LO a temporal profile in 
which matches that of the vacuum modes (as can be seen from the correlation functions). In 
this case the normalized power measured by a Fourier-limited spectrum analyzer integrating 
over a time large with respect to 1 /70 is given by 



/"/J(i J . rt+-lo rt+T 
7T & / dT'e~^-^E L o{T)E L o(T')C'i(ry) 

-7r/T 27T Jt Jt 
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again with a symmetrical expression for the second homodyne detection. M and S are integrals 
depending on To, 70 and C . M represents the shot noise (i.e. the noise level for uncorrelated 
atoms, and, therefore, uncorrelated fields), so that S/N is the signal-to-noise ratio of the 
measurement. M can be easily measured in a preliminary experiment in which the atoms are 
prepared in a coherent spin state. The signal-to-noise ratio can be shown to be close to 1 
when 70 To ^> 1 and C ^S> 1 [15]. Note that this is possible because we have chosen the right 
matching profile for the LO. For short times the atomic entanglement is given by 



^[Pi(0) + P 2 (0)] 



^[A 2 (J,i 



J x2 ) + A 2 {J yl + J y2 )} =l a 
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We indeed retrieve the value of the inseparability criterion JfJJl for the atomic ensembles, i.e. 
measuring the atomic entanglement is equivalent to measuring the entanglement between the 
outgoing modes. It is then possible with this technique to detect the atomic entanglement 
with nearly 100% efficiency. 

In conclusion we have proposed a scheme to achieve continuous entanglement between 
atomic ensembles in cavities using a pair of EPR-correlated beams. The entanglement can 
be stored for a long time in the ground-state atomic spins and retrieved at will in the fields 
exiting the cavities by switching on the control fields. We also propose a technique to perform 
the atomic entanglement readout in a single shot measurement with one local oscillator. It is 
worth noticing that all the results obtained in an EIT configuration can be readily transposed 
in a Raman configuration, in which a perfect entanglement storage is also predicted in the 
regime 70 <C (1 + 2C)Tr <C k, 7 (where Tr = jtt 2 /A 2 is the Raman optical pumping rate). 
We would like to point out that, although the storage time is given by the inverse of the natural 
decay rate of the ground state I/70, the memory bandwidth, i.e. the frequency bandwidth over 
which the entanglement is stored is much broader (oc 70), because of the cavity interaction 
with the field. Such a quantum memory could allow to store and manipulate entanglement 
which is a key challenge for quantum communication and information. 

* * * 
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